Search for: All records

The early Martian atmosphere had 0.25 to 4 bar of CO₂but thinned rapidly around 3.5 billion years ago. The fate of that carbon remains poorly constrained. The hydrothermal alteration of ultramafic rocks, rich in Fe(II) and Mg, forms both abiotic methane, serpentine, and high-surface-area smectite clays. Given the abundance of ultramafic rocks and smectite in the Martian upper crust and the growing evidence of organic carbon in Martian sedimentary rocks, we quantify the effects of ultramafic alteration on the carbon cycle of early Mars. We calculate the capacity of Noachian-age clays to store organic carbon. Up to 1.7 bar of CO₂can plausibly be adsorbed on clay surfaces. Coupling abiotic methanogenesis with best estimates of Mars’ δ¹³C history predicts a reservoir of 0.6 to 1.3 bar of CO₂equivalent. Such a reservoir could be used as an energy source for long-term missions. Our results further illustrate the control of water-rock reactions on the atmospheric evolution of planets. 

The continental lithospheric mantle plays an essential role in stabilizing continents over long geological time scales. Quantifying spatial variations in thermal and compositional properties of the mantle lithosphere is crucial to understanding its formation and its impact on continental stability; however, our understanding of these variations remains limited. Here we apply the Whole-rock Interpretive Seismic Toolbox For Ultramafic Lithologies (WISTFUL) to estimate thermal, compositional, and density variations in the continental mantle beneath the contiguous United States from MITPS_20, a joint body and surface wave tomographic inversion for Vp and Vs with high resolution in the shallow mantle (60–100 km). Our analysis shows lateral variations in temperature beneath the continental United States of up to 800–900 °C at 60, 80, and 100 km depth. East of the Rocky Mountains, the mantle lithosphere is generally cold (350–850 °C at 60 km), with higher temperatures (up to 1000 °C at 60 km) along the Atlantic coastal margin. By contrast, the mantle lithosphere west of the Rocky Mountains is hot (typically >1000 °C at 60 km, >1200 °C at 80–100 km), with the highest temperatures beneath Holocene volcanoes. In agreement with previous work, we find that the chemical depletion predicted by WISTFUL does not fully offset the density difference due to temperature. Extending our results using Rayleigh-Taylor instability analysis, implies the lithosphere below the United States could be undergoing oscillatory convection, in which cooling, densification, and sinking of a chemically buoyant layer alternates with reheating and rising of that layer. 

Abstract To quantitatively convert upper mantle seismic wave speeds measured into temperature, density, composition, and corresponding and uncertainty, we introduce theWhole‐rockInterpretativeSeismicToolboxForUltramaficLithologies (WISTFUL). WISTFUL is underpinned by a database of 4,485 ultramafic whole‐rock compositions, their calculated mineral modes, elastic moduli, and seismic wave speeds over a range of pressure (P) and temperature (T) (P = 0.5–6 GPa,T = 200–1,600°C) using the Gibbs free energy minimization routine Perple_X. These data are interpreted with a toolbox of MATLAB® functions, scripts, and three general user interfaces:WISTFUL_relations, which plots relationships between calculated parameters and/or composition;WISTFUL_geotherms, which calculates seismic wave speeds along geotherms; andWISTFUL_inversion, which inverts seismic wave speeds for best‐fit temperature, composition, and density. To evaluate our methodology and quantify the forward calculation error, we estimate two dominant sources of uncertainty: (a) the predicted mineral modes and compositions, and (b) the elastic properties and mixing equations. To constrain the first source of uncertainty, we compiled 122 well‐studied ultramafic xenoliths with known whole‐rock compositions, mineral modes, and estimatedP‐Tconditions. We compared the observed mineral modes with modes predicted using five different thermodynamic solid solution models. The Holland et al. (2018,https://doi.org/10.1093/petrology/egy048) solution models best reproduce phase assemblages (∼12 vol. % phase root‐mean‐square error [RMSE]) and estimated wave speeds. To assess the second source of uncertainty, we compared wave speed measurements of 40 ultramafic rocks with calculated wave speeds, finding excellent agreement (V_pRMSE = 0.11 km/s). WISTFUL easily analyzes seismic datasets, integrates into modeling, and acts as an educational tool. 

Abstract The India‐Eurasia collision is a key case study for understanding the influence of plate tectonic processes on Earth's crust, atmosphere, hydrosphere, and biosphere. However, the timing of the final India‐Eurasia continental collision is debated due to significant uncertainty in the age of the collision between the Kohistan‐Ladakh arc (KLA) and Eurasia along the Shyok suture zone. Here we present paleomagnetic results that constrain the Karakoram terrane in northwest India to a paleolatitude of 19.9 ± 8.9°N between 93 and 75 million years ago (Ma). Our results show that the Karakoram terrane was situated on the southern margin of Eurasia in the Late‐Cretaceous. Our results indicate that the KLA and Eurasian continent had a not converged until <61.6 Ma, placing a Paleocene older limit on the age of final closure of the Shyok suture zone. This suggests that the India‐Eurasia collision in northwestern India likely occurred after the closure of the oceanic basin between the KLA and Eurasia. The Paleocene collision event affecting India that has been widely interpreted to represent final India‐Eurasia collision instead records the arc‐continent collision between the KLA and the northern edge of India prior to final India‐Eurasia collision. Final India‐Eurasia collision in northwest India most likely occurred after the closure of the oceanic basin between the KLA and Eurasia.